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iJS ' Abstract 

> : 

^^ , We study Higgs boson production from gluon-gluon fusion at mid-rapidity in high energy proton- 

ON ■ nucleus colhsions.For this process the presently still little known gluon distribution function h^ ^{x, k±) 

^N I might give a numerically relevant contribution. We show by explicite calculation that using CGC 

i__i' (color glass condensate) model input the result obtained in the naive kt factorization approach 

f~i I matches the result obtained in the TMD factorization framework for a dilute medium. We also ver- 

Mh' ify the earlier finding [14] that the k^ factorization formalism for Higgs production breaks down in 

Oj! a dense medium. In doing so we formulate a hybrid model which allows one to treat such reactions 

D ■ theoretically. 

^ ■ 1 Introduction 

>: 

'nI" ' In recent years the theoretical understanding of transverse momentum dependent (TMD) parton 

distributions has made tremendous progress, up to the point that their further investigation can serve 

as one of the motivations for a new electron-ion collider [1]. Still there are many aspects which need 

^r^ I further study. An especially interesting one is relevant for Higgs-production at the LHC. In this 

^D ■ case it was realized that for pp collisions in addition to the standard gluon contribution one gets a 

^_j , contribution proportional to a little known distribution function, usually refereed to as the distribution 

;^ I of linearly polarized gluons [2,3] {h^^ in the notation of Ref. [4]). Therefore, h^^ has attracted a lot 

.^H ' of attention recently. This new distribution function is the only spin dependent gluon TMD for an 

/\ . unpolarized nucleon/nucleus, and may be considered as the counterpart of the quark Boer-Mulders 

C^ ! function hi'^{x,k±) [5]. However, in contrast to the latter, h^^ is time-reversal even implying that 

initial/final state interactions are not needed for its existence [6,7]. This distribution function is of 

phenomenological interest, especially for small- a; physics at RHIC and LHC because a calculation in 

the saturation model [8] showed that its contributions are (at small-x) as large as those proportional to 

the unpolarized gluon distribution. Fortunately, it has been shown that h^ ^ can be accessed, at least 

in principle, through measuring, e.g., azimuthal cos 20 asymmetries in processes such as jet or heavy 

quark pair production in electron-nucleon scattering as well as nucleon-nucleon scattering. Other 

promising observables are cos2i?i) asymmetries in photon pair production in hadron collisions [9-11]. 

Such measurements should be feasible at RHIC, the LHC, and a potential future Electron Ion Collider 

(EIC) [1,12] and could play an important role to establish saturation effects. More recently, it has been 

found that the linearly polarized gluon distribution may affect the transverse momentum distribution 

of Higgs bosons produced from gluon fusion for pfj± ^ mn, where ph± and mn are the Higgs 

transverse momentum and mass respectively [13, 14]. The authors of Ref. [13] proposed that the 

effect of linearly polarized gluons on the Higgs transverse momentum distribution can even be used, in 

principle, to determine the parity of the Higgs boson experimentally. Transverse momentum dependent 

factorization has been re-examined by taking into account the perturbative gluon-radiation correction 






to h^^ [14]. The complete TMD factorization results for Higgs boson production are consistent with 
earlier findings based on the Collins-Soper-Sterman (CSS) formalism [17] and soft-coUinear-effective 
theory [18]. Also, the transverse momentum resummation formalism applied to di-photon production 
in pp collisions [19] is closely related. 

Besides their obvious phenomenological interest, these investigations of Higgs production are also 
interesting from a more theoretical point of view. The theoretical description of transverse momentum 
dependent processes unavoidably includes gauge links in one or the other form. While the starting 
expressions for the different approaches look often quite different with respect to these gauge links, 
it turned out that quite often the resulting cross sections can be mapped onto one another in some 
approximation and for a suitable kinematic window. In the following we will discuss three such ap- 
proaches, TMD factorization, kt factorization and a new hybrid approach. 

TMD factorization 

We will use the term TMD factorization in the sense of [20], which defines hard and soft factors, such 
that, e.g., the Higgs boson production cross section with Pi <C M reads [14] 

X \^xig{xi,ki^±) X2g{x2,k2,±) S{£±,fip) H{M'^,^ip) 

+ 1,2 1.2^ 1 Xih\;^{xi,ki^i_) X2h'\;^{x2,k2,L) 



X Sh{^±,^lp)Hh{M\|lp)] 



[I) 



with the soft factors S and Sh and the hard interaction factors H and Hh. TMD factorization is 
probably the formally most complete and reliable scheme, but often also the calculational most de- 
manding. For specific questions other schemes might be more economic. For example for all-order 
proofs of factorization SCET is a promising alternative [15] and for qualitative phenomenological anal- 
ysis TMD factorization promisees a substantial simplification. Within TMD factorization there also 
exist different approaches. To be specific we use this term for the formulation of Collins et al. for 
which it is crucial to define the gauge links slightly off the light-cone. Within SCET it is possible to 
keep the gauge links on the light-cone [16]. In principle both approaches should give consistent results 
for physical observables when expanded in an appropriate manner, but it is non-trivial to map, e.g., 
evolution in both schemes onto one another. 

kt factorization 

In contrast, the naive kt factorization scheme invokes some approximations. In this formulation there 
is no linearly polarized gluon distribution function which is equivalent to the statement that it has 
to have the same functional form as the normal unpolarized gluon distribution such that it cannot 
be discerned. This fact demonstrates clearly that kt factorization can in general not be a good 
approximation because within CGC framework both gluon distributions could differ substantially for 
Qs ^ k±_ as was first noticed in [8]. This can be best seen from the following expressions for the 
Weizsacker- Williams (WW) unpolarized gluon distribution denoted by Gww{x^k±) and the WW 
type linearly polarized gluon distribution /i^ f^/^^(x, A;j_) derived in the CGC formalism [8,24,25], 



.G,„„..,J = ^_^/,V,e-.S.-.-^(l 



±9 ,^^ -^ ^ ^c - 1 o /■ .„„ U^l 



Here S± is the transverse area of the target nucleus, and k± = \k±\. Q^ = agNcf^A^T^ ^/(^I^qcd 
is the gluon saturation scale with fiA being a common CGC parameter. Note that our convention for 
h^ ^r\Y differs from that in Ref. [8] by a factor 1/2. In general, both gluon distributions are different 
though they become identical in the dilute region, i.e. for k_i » Qg 

hybrid approach 

Our main strategy is to perform the calculations for proton-nucleus reactions in an approach where 
the nucleus is treated in the CGC framework [28-31], which effectively takes soft gluons into account, 
and the proton in the so-called Lipatov approximation [32-35] . We restrict ourself to studying Higgs 
production in the plain MV model [28] without considering small x evolution effects although the 
latter could be done in principle [36-42] by solving the general JIMWLK evolution equation [43] for 
quadrupole operator. Neglecting evolution should be a good approximation because the MV model is 
valid in the range x « 0.01 — 0.1, which is the relevant kinematical regime for Higgs boson production 
at LHC, while radiative corrections only become important below a certain scale x ~ 0.01. 

This article is organized as follows. In Sec. II, we introduce the hybrid approach and use it to 
reproduce the well known result for soft gluon production in pA collisions. The Higgs boson production 
in pA collisions is computed in the same approach in Sec HI. We also demonstrate consistency between 
the results obtained in TMD factorization and kt factorization in the dilute region within the CGC 
model. We summary our paper in Sec. IV. 

2 Soft gluon production in proton-nucleus collisions 

In this section, we introduce a hybrid approach and reproduce the well known result for gluon pro- 
duction in high energy proton- nucleus scattering [44-49]. In the next section we shall generalize this 
approach to Higgs production. Let us first consider the general case of soft gluon production 

A{Pa)+p{Pb)^9{1)+X . (4) 

We assume that the nucleus is moving with a velocity very close to the speed of light into the positive 
z direction, while the proton is moving in the opposite direction. It is convenient to use light-cone 
coordinates for which P^ = P'^p^ and Pg = P^n^ with p = (1,0,0,0) and n = (0,1,0,0). The 
corresponding partonic subprocess is represented by g A{ki) + gp{k2) -^ g{l), where k^ = xiP^p^ + k^^^ 
denotes the total momentum carried by multiple gluons from the nucleus, and kt^ = X2P^n^ + kl^^ is 
the momentum of the gluon from the proton. We chose to work in the light-cone gauge of the proton 
{A~ = 0), for which the polarization tensor of a produced gluon is given by, 

J2 e^e*' = -g^'' + ^ ^^ . (5) 

As mentioned above, to facilitate our calculation, a hybrid strategy has been adopted, in which the 
nucleus is treated in the CGC model, while on the side of the dilute projectile proton one makes 
the so-called Lipatov approximation [32-35]. At small x the gluon radiation cascade shows a strong 
ordering in rapidity. Or in other words, color source carries much larger rapidity than that radiated 
gluon does. It has been shown that a fast moving color source can be treated as eikonal line in the 
strongly rapidity ordered region. The operation of introducing these eikonal lines is refereed to as 



k 

Figure 1: Feynman rules for the eikonal line, which is represented by a double line, a, b, c denote color 
indices. 

Lipatov approximation [35]. Its validity has been confirmed also by solving the classical Yang- Mills 
equation [47-49]. In these calculations, the gluon field induced inside a proton by a weak color source 
and weak color source itself are treated as a small parameter when solving classical Yang-Mills equation 
perturbatively. An analytic solution for the gauge field was obtained in lowest order of the incoming 
gluon field in various gauges [47-49]. 

For the process of gluon production in pA collisions, the relevant eikonal line is the past-pointing 
one which is built up through initial state interactions between the color source inside the proton and 
the background gluon field. The interaction between the classical gluon field and the final state gluon 
emitted from the color source inside the proton does not change this general statement because the 
imaginary part of the scattering amplitude cancels between the different cut diagrams once the final 
states are integrated out. The prescription to treat the eikonal propagator is fixed by this choice. The 
relevant Feynman rules, illustrated in Fig. 1, were given in Ref. [35]. Note that the prescription for 
past-pointing eikonal propagators differs from that for future-pointing eikonal lines. 

It is worthwhile to mention that to preserve gauge invariance one has to take both, gluon fusion 
and the interaction between the color source inside the proton and the strong classical gluon field of 
the large nucleus into account. This is due to the fact that the incoming gluon from the proton is 
off-shell and off-shell quantities are, in general, not gauge invariant. Both types of interaction are 
shown in Fig. 2. 

The multiple scattering between incoming gluon (or eikonal line) and the classical color field of 
the nucleus can be readily resumed to all orders as has been done for the scattering of a quark by a 
background gluon field [51,52]. The resumed multiple scattering gives rise to a path-ordered gauge 
factor along the straight line that extends in x~ from minus infinity to plus infinity. More precisely, 
for a gluon or eikonal line with incoming momentum being k and outgoing momentum being k + q, 
the path-ordered gauge factor reads, 

27r%-K[t/-l](Qx), (6) 

with 

[U - l]{q±) = j d^x^e-'^i^-^^pix^) - 1] , (7) 

and 

U{x^) = (7'e-''^/-r'^--^+(--' -^))^ , (8) 

where A'^ = A^t'^ is the gluon potential in the adjoint representation and t^j^ = —ifabc ■ 

We use this as a building block to compute the amplitude for gluon production in high energy pA 
collisions. It is easy to verify that the contribution from the third diagram vanishes because both k[ 
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(6) 




Figure 2: The diagrams contributing to gluon production. The dash hues represent the re-summed 
interactions of the incoming gluon or color source inside the proton with the classical color field of the 
nucleus. 

poles are located in the same half plane. Consequently, we are left with the contributions from the 
first two diagrams. The calculation of these two diagrams is straightforward. Collecting all pieces, the 
differential cross section reads, 



da 
d?ldy 



IT 



2d^ki± 



{N^-l)P J (2^)3^^?±^2<7(:r2,A:2±)[C/-l](A:i±)[C/^-l](A;ix) 



d^ki±_ 

(iV2 - l)/i J 7^ 



kl^X2g{x2, k2±)U{ki±)U\ki±) 



(9) 



where g{x2,k2± = |^2±| = |^± — ^i±l) denotes the un-integrated gluon distribution of a proton, and 
y is the rapidity of the produced gluon. The factor 2/(27r)'^ associated with phase space integration 

is chosen such that for single gluon target, / (27rV^'^ ^^(^(^ix)?^ (fei±))giuon = xi5{l — xi) at lowest 
non-trivial order (see, for example, Ref. [30]). To obtain the above result, we have defined the normal- 
ization factor and the flux factor to be fc2_|_/(2A:2 ■ p{N^ — 1)) and l/(2/c2 -p), respectively, rather than 
kfj_k2j_/ {4:XiX2Pa ■ Pb{^c ~ 1)^)) ^/{^xiX2Pa • Pb) used in Ref. [35], since the Lipatov approximation 
is only justified for the proton side. 

The next step is to compute the expectation value of a Wilson line in the plain McLerran- 
Venugopalan model [28]. By averaging over color sources with a Gaussian distribution, one finds [53- 
55], 

(TrC/(fix+rx)C/t(ii^)\ =(iV2_i)exp|^^| . (10) 



To proceed further, one notices that the dipole type gluon distribution in the adjoint representation 
is given by [46,56,57], 



xGdp{x, k^) = S^A:2 / ^ e-'^^-^ e ^ 



27r2a, ^J (27r)2 

With the help of the above two equations, the differential cross section can be written as, 

da An'^a.^Nr , ,^, , , n ^ / , x 

d ki_i_X2g{x2,k2±)xiGDp{xi,ki±) 



dHdy (Af2 _ 1)/ 



± 



(11) 



(12) 



which is in full agreement with the results of earlier calculations [44-50]. To see this, one has to iden- 
tify nx2g{x2,k2±) = d[x2g{x2,Q'^)]/dQ'^\Q2^i,2^ and irxiGDp{xi,ki±) = d[xiG{xi,Q'^)]/dQ'^\Q2^q^, 



where g{x2,Q'^) and G{xi,Q'^) are the integrated gluon distributions of proton and nucleus. In par- 
ticular, Eq.[12] demonstrates that the gluon production cross section can be expressed in terms of 
the gluon distribution in a rather straightforward manner. In the next section, we apply this hybrid 
approach to calculate the Higgs boson production in pA collisions. In contrast to gluon production for 
which the dipole gluon distribution appears in the cross section, one has to use Weizsacker- Williams 
gluon distributions in the calculation for Higgs boson production, as we show below. 

3 Higgs boson production in proton-nucleus collisions 

Now we turn to Higgs boson production in proton-nucleus collisions. Let us start by introducing the 
matrix element definition for gluon TMDs that generate the Higgs transverse momentum distribution 

[2,4,58,59], 

= -YXiGww{xi,ki±) + (k\±klj_ - -6f]xlhfl^r^^{xl,kl±) , (13) 

where k\j_ = k\j_/ki±, and 6^ = —g^-' + {p^n^ + p'n'^)/p ■ n. The gauge link extends to the past: 

Ly = V e~^^ J oo- ^ ^^ '^^^ The two leading power gluon TMDs ^^^/^^(xi, A;i^), /ij^ ^jy(xi, /ci^) are 
the usual WW type unpolarized TMD gluon distribution and WW type TMD distribution of linearly 
polarized gluons respectively. As we show below, both gluon TMDs contribute to the differential cross 
section for Higgs production. 

Higgs boson production through gluon-gluon fusion has been studied in the context of the kf 
factorization formalism [21-23]. However, the authors of Ref. [14] argued that kt factorization can 
break down in a dense medium where the WW type linearly polarized gluon distribution is different 
from the usual gluon distribution. Then the Higgs production cross section cannot be expressed only in 
terms of the usual gluon distribution and TMD and k± factorization give different results. However, 
they also argue that one should be able to modify kt factorization to establish an effective TMD 
factorization at small x. One possibility to do so was recently proposed in Refs. [26,27]. Generally 
speaking, at very small x higher twist contributions are as important as the leading twist ones because 
of the high gluon density. Therefore, in order to arrive at the mentioned effective TMD factorization, 
an analysis including all higher twist contributions is crucial. For the unpolarized case it was shown 
that the results for two-particle correlations in high energy scattering using the proposed effective 
TMD factorization are in agreement with the results obtained by extrapolating the CGC calculation 
to the correlation limit [26, 27], where the transverse momentum imbalance between the two final 
state particles (or jets) is much smaller than the individual transverse momenta. By applying a 
corresponding power counting in momentum space in the correlation limit, a complete matching 
between the effective TMD factorization and the CGC formalism has also been found in the polarized 
case [8]. Later, a calculation performed in position space led to the same results [27]. 

Inspired by Ref. [14], we carry out an explicit calculation for Higgs boson production in pA collisions 
using the CGC formalism and verify the conjecture that the effective TMD factorization and CGC 
approach provide the same result in the dense medium region, while the kt factorization is only valid 
in the dilute region. The starting point is the effective Lagrangian for Higgs boson production, 

Ceff = -\g^^F^,,F'^^'^ , (14) 
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Figure 3: The generic diagrams contributing to Higgs boson production. The multiple scattering of 
the incoming gluon or color source inside the proton by the classical gauge field have to be re-summed 
to all orders, pn denotes the Higgs boson momentum. 



which is valid in the heavy top quark limit, where $ is the scalar field and F^'^ the gluon field 
strength, g^p is the effective coupling. The same effective Lagrangian has also been used to study 
gluon saturation in semi-inclusive DIS off large nuclei [45]. From the above Lagrangian, we can read 
off the basic vertices for the Higgs boson coupling to gluon. The corresponding Feynman rule for a 
Higgs boson coupling to two off-shell gluons carrying the momenta ki, k2, and color indices a, b is 
given by 



ig^^ab (g^^h ■ k2 - k'^k 



(15) 



We argue that Higgs boson production through a multiple gluons fusion process is not enhanced by 
saturation effects for the following reasons: First, the dominant component of the classical gauge 
field A^ is generated by the color sources inside the nucleus in a reasonably local way [45]. Second, 
the coherence length for Higgs boson production is very small due to the large top quark mass. 
Therefore, to calculate the transverse momentum dependent cross section for Higgs boson production 
it is sufficient (in the saturation region) to use the effective vertex given above. 

The relevant diagrams for Higgs boson production are shown in Fig. 3. Before the incoming gluon 
from the proton combines with a gluon from the classical gauge field of the nucleus to form a Higgs 
boson, this gluon has initial state interactions with background gluons as it passes through the nucleus. 
Although they do not modify the integrated production rate, the initial state interactions change the 
transverse momentum distribution of the produced Higgs boson. The initial state interaction between 
the color source inside the proton and the classical gluon field of the nucleus should in principle also 
be taken into account. However, such initial state interaction, shown in Fig. 3(c) and Fig. 3(d) vanish 
because the k'l poles are in the same half plane. Only Fig. 3(a) and Fig. 3(b) give non-vanishing 
contributions. Resuming gluon re-scattering to all orders, as illustrated in Fig. 3(b), and combining 
it with the contribution from Fig. 3(a), one obtains the production amplitude, 



M =g^{k2-p) 



fk[^fk2±-k\ 

(2vr)2 



i± 



A-2 



{L{ki^ - k[^)A+{k[^))AM^2,k: 



'2±) 



(16) 



where the ki' and k^ components have been integrated out. k2 = ph — ki denotes the momentum of 
the gluon from the proton with pn being the Higgs momentum. (j)p{x2, k2±_) represents the probability 
amplitude for finding a gluon carrying a certain momentum inside the proton, with X2g{x2,k2±) = 
4>p(t)p- It is convenient to introduce the gauge potential in impact parameter space. 



A+{k[ 



i±) 



dy d y±e' 



nP+y 



ik' 



^■y^A+{y~,y^) 



(17) 



Replacing the exponential e*^i ^ by e*^^ ^ in the above formula is justified in the leading loga- 
rithm approximation. Moreover, one notices that L{ki_i — k[j_) is a Wilson line which starts from ^~ 
being minus infinity and ends at the space-time point y~ , 



We proceed by partial integration and by performing the integral over /c^_|_. The amplitude then can 
be written as, 






Using the same normalization and flux factors as in the previous section, the differential cross section 
becomes, 

d^a TTgl f2d'^ki±k2±,i^2±,j , , ^ /" ,3 ,3 

- ' -X2g{x2,k2±) d-'rd-'y 



cPpH^dy 32 J (27r)3 kl 



2± 



X ^^^^P^^^-+^>^-^-^^ {F+\r- + y-,r^+ y^) L^y Ly F^^{y',y^))A 



= ao d^ki±X2g{x2,k2±)xi Gww{xi,ki±) + (2{ki_i_ ■ k2±f - Ij hf^ly^{xi,ki_i_) (20) 

where (Tq = -g^ is the leading-order cross section for scalar-particle production from two gluons. Here, 
y and Ph± are rapidity and transverse momentum of the Higgs particle. In the second step of above 
derivation, we have made use of the normalization conditions for the average over the CGC wave 
function: (1)^ = 1; and for the nuclear state \P) carrying momentum P: {P'\P) = 2P^{2tt)'^6{P'^ — 
P'+)6'^{P^-P'^) [27]. 

As observed in Ref. [14], one automatically takes into account the contribution from the linearly 
polarized gluon TMD in kt factorization. In other words, the usual unpolarized gluon distribution 
of the proton is the same as its linearly polarized gluon distribution in the Lipatov approximation. 
By noticing this fact, one finds that the differential cross section computed in TMD factorization [13, 
14] completely agrees with that derived in the CGC approach. It is also worthwhile to mention 
that the gluon distributions entering the cross section for Higgs boson production are the WW type 
distributions as expected. Furthermore, in order to compare this with results from the kt factorization 
formalism, we recall that the WW type distribution Gww and h-^ y^y^ become identical in the dilute 
region where ki± ^ Qs, but differ in a dense medium [8]. Therefore, in the dilute region, the 
differential cross section can be simplified to 

do / d^ki^X2g{x2,k2i_)xiGww{xi,ki^)2{ki^ ■ k2±f (21) 



(PpH±dy 



which agrees with the well known result obtained from the kt factorization approach [21, 22] at low 
Higgs transverse momentum. Thus we conclude that also for Higgs production in pA the kt factoriza- 
tion formula is only valid in the dilute region. Similar conclusions have also been drawn for r/' meson 
production [60] and heavy quark pair production [54] in pA collisions. 



4 Suminary 

We developed a hybrid approach for calculating particles production at central rapidity in pA collisions, 
in which the dense target nucleus is treated in the color glass condensate model, while on the side of 
the dilute projectile proton the Lipatov approximation was used. As a test of the method, we first 
reproduced the well known result for soft gluon production in pA collisions using this hybrid approach. 
Then, we derived the differential cross section for Higgs boson production from gluon fusion in pA 
collisions. It turned out that the result obtained in our hybrid approach is completely equivalent to 
that computed in TMD factorization. In the low-density limit of pA collisions, we also recover the 
result of the naive kt factorization formalism that describes Higgs boson production in pp collisions 
adequately. 

The approach developed in this article also can be applied to study the production of other color- 
neutral particles or heavy quark pair production in pA collisions. One may expect that the CGC 
formalism and TMD factorization will yield the same results for these processes in a certain kine- 
matical region. As a consequence, the Weizsacker- Williams and dipole type linearly polarized gluon 
distributions could be extracted by measuring azimuthal asymmetries in these processes. We will 
address these issues in a forthcoming paper [61]. 
Acknowledgments: This work has been supported by BMBF (OR 06RY9191). 
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